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Abstract 

The 5-matrix of the Ising Model can be obtained as particular limit of the 
roaming trajectories associated to of the S-matrix of the Sinh-Gordon model. Using 
the form factors of the Sinh-Gordon, we analyse the correspondence between the 
operators of the two theories. 



1 Introduction 



Given an elastic factorized S-matrix of a 2-D system with a mass scale M, we can cal- 
culate its ground state energy Eq(R) = —irc(MR)/6R on an infinite strip of width R, 
by means of the Thermodynamical Bethe Ansatz (TBA) fl], ^ . For those models where 
the S'-matrix has a well-defined field theory correspondence || £|, §], the scaling func- 
tion c(MR) has a smooth behaviour, monotonically decreasing from the limit value c(0) 
(where it coincides with the effective central charge of the CFT of the ultraviolet limit) 
to c(oo) = (which corresponds to massive regime). However, since the TBA only em- 
ploys an S'-matrix without questioning its field theory interpretation, it can be also used 
to investigate the finite-size behaviour of any quantum theory axiomatically defined in 
terms of a scattering amplitude, provided it satisfies the usual constraints of unitarity and 
crossing symmetry. From this point of view, ALB. Zamolodchikov proposed in ref. || a 
simple purely elastic scattering theory which under the TBA analysis reveals a remarkable 
finite-size behavior. Such theory contains a single particle bosonic state with mass M and 
two-particle scattering amplitude given by 

smh^j-|cosh^ 
[P) sinh/3 + icoshA, ' 1 J 

where (3q is a real parameter. S(/5) has two simple poles in the unphysical sheet at 
positions (3 = — y ± /3 which correspond to a resonance particle. The presence of a 
scale (3q for real values of the rapidities drastically influences the finite-size behaviour of 
the model. In fact, solving numerically the TBA equations associated to the S-matrix 
(|TTT|), for sufficient large values of j3o, c(r) develops a "staircase" pattern with a series 
of plateaux at the discrete values c = 1 — G° = 3,4, . . .) which coincide with the 

central charges of unitary minimal models M. p [0, ||. Hence the Roaming Trajectory 
Model (RTM) is a one-parameter family of Renormalization Group flows interpolating 
between all the fixed points M. v : each trajectory starts from the limiting fixed point 
M.OO and then proceeds on the critical surface through the hopping M. v — > M. v -\ until 
it arrives in the neighborhood of the fixed point After this last step, it develops a 
finite correlation length and gives rise to the usual massive infrared behaviour. From the 
TBA analysis it also follows that the roaming trajectories spend approximately the same 
fraction (3q of the Renormalization Group "time" x = log MR/2 near each fixed point, 
therefore making more pronounced the multiple crossover phenomena for large values of 
(3q. Although a local interpretation of the RTM has been pursued in terms of conformal 
perturbation of the models M. v visiting along the trajectories |TD[, it is worth to consider 
the RTM as special analytic continuation of the Sinh-Gordon model in such a way to take 
advantage of the recent exact solution of this model [|ll|, Purpose of this letter is 
to show, as simplest application of this idea, how to relate the operator content of the 
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Sinh-Gordon model to that of the Ising model which is the first jump in the staircase 
series. 



2 The Sinh-Gordon model 



2.1 Main features 



The Sinh-Gordon Model (SGM) is a 2-D Affine Toda Field Theories [0 with one bosonic 
field <f)(x) and bare action given by 



A = d 2 x 



1~ ^ Ml 



JPu$) r cosh^0(x) 

^ 9 



(2.1) 



The integrability of the model permits the determination of the factorizable elastic S- 
matrix which is given by 



sinh 3 — % sin ^ 
S(P,B)= , ^ , ^ , 2.2 
sinh p + 2 sin ^ 



where -B(^) = 8 ^ 2 • F° r rea l values of the coupling constant g, the S"-matrix has no poles 
in the physical sheet and consequently no bound states, but on the contrary it presents 
two zeroes at the crossing symmetric positions inB/2 and Z7r(2 — B)/2. It is easy to see 
that in the analytical continuation B — > 1 ± — 3q the zeros move along a direction parallel 
to the real /3-axis and the 5-matrix (|2.2| ) exactly coincides with the scattering amplitude 
of the RTM ||]. This observation becomes useful in the light of the fact that the SGM 
has been recently solved by computing the matrix elements of local operators. 

2.2 Form Factors 

A complete knowledge of a QFT is encoded into the matrix elements of local operators 
Ok on the asymptotic states, the so-called Form Factors (FF) |15| 



F n fc (A,...,/3 n ) =<0\O k (0)\p 1 ,...,p n > . (2.3) 
In the case of the SGM at real coupling constant, the FF of local scalar operators have 



been determined in [11, 12 1. We briefly recall their main properties, referring the reader 



to the original references for their detailed discussion. They can be parameterized as 

i<j l x « I" x j) 
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where Xj = e ft and (3^ = A — /%. F m i n (/3) is an analytic function given by 
F min ((3,B) = Af(B) E(J3,B) 
Z{P,B) = exp 



rf x sinh(f ) sinh (f (1 - f )) sinhf 



o x 



sinh x 



sin 



x£ 
2tt 



(2.5) 



JV(-B) = exp 



d x sinh (2^) sinh (|(1 



2 / y sinh f 



X 



sinh x 



(P = in — (3). F min (P,B) has a simple zero at the threshold P = and no poles in the 
physical strip < Im/3 < ir, with an asymptotic behaviour lim^oo F min (P, B) = 1. In 
eq. ( |2.4| ) if ^ are normalization constants which depend on the operator one is considering. 
The functions Q„( symmetric polynomials in the variables Xj, solutions of 

the recursion equations which link the n-particle and the (n + 2)-particle form factors 



l lim (/? - /?)F„ fc +2 (/? + iTT, /?, Pi, p 2 ,...,Pn 

13^/3 



(2.6) 

For FF of spinless operators, the total degree of Q k n is equal to n(n — l)/2 whereas their 
partial degree in each variable X, depends on the operator Ok which is considered. It 
was shown in ref. [|12| that a general solution for the Q\ can be written in terms of the 
so-called elementary solutions Q n {p) given by[| 



Q n (p) = detiWy(p) 



(2.7) 



where Mij(p) is an (n — 1) x (n — 1) matrix with entries Mij(p) = o~2i~j [i — j + p] {o~i are 
the elementary symmetric polynomials [ITH] and p an arbitrary integer). 



2.2.1 Form Factors of 0(x) and 9(x) 

Important operators of the SGM are the elementary field <p(x) and the trace of the stress- 
energy tensor O(x). They are odd and even operators respectively under the Z 2 symmetry 
of the model with normalizations given by < | 0(0) | P >= 1 and < P | O(0) | P >= 
2ttM 2 , where M is the physical mass. The whole set of FF of the elementary field 4>(x) 
is given by 

They are automatically zero for even n (in agreement with the Z 2 parity of the model) 
whereas for odd n they vanish asymptotically when — > oo, as follows from the LSZ 

1 We have suppressed the dependence of Q(p) from the variables Xi. 
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reduction formula. Concerning the FF of Q(x), -Fjn+i = whereas Ff n are given by 



F 2 „(/5 l5 . . . ,/? 2n ) - ^ ^ ) Q 2 n(l) 11^—7- , (2-9) 

and they go to a constant when $ — > 00 
2.2.2 Kernel Solutions 

The general structure of the FF of the SGM is that of a sequence of finite linear spaces 
whose dimensions grow linearly as n increasing the number 2n — 1 or 2n of external parti- 
cles. In fact, at each level of the recursive process the space of the FF is enlarged by includ- 
ing the kernel solutions of the recursive equation (|2.6| ), i.e. Q n (—x,x,Xi, . . . , av- 2 ) = 0. 
With the constraint that the total order of the polynomials is n ^ n ~ 1 ^ ; the kernel is unique 
and given by Ti n (xi, . . . , x n ) = det cr 2 j_j. This solution gives rise to the half-infinite chain 
under the recursive pinching x\ — —X2 = x 

• • • - Q { & - Qui - Q { : } = - (2.10) 

and therefore the whole space of FF presents the foliation structure^] 

• • • ~~ * Qn+4 ~ > Qn+2 ~~ > QtP ~~ * Qn-2 ~^ ■ ■ • ~~ > Q?.^ 



!n+2 ~~* ^in Vn-2 • • • ~^ V3 

Qn+4, ~~ * Qn+2 ~~ * Qn^ ~^ Qn- 



Vn+4 ' ' 



(2.11) 







(n+2) v 
n+4 ~~ * ^n+2 



The explicit expressions of such solutions can be found by determining the linear combi- 
nation of Q n (k) which reduces to S n at the level n. 



3 Violation of the c-theorem sum rule in the RTM 

Since the RTM may be seen as the SGM at B = 1 ± — /3 , it is natural to study the 
behaviour of the FF of the latter model under this analytic continuation. As we show, 
the presence of a scale (3q in the rapidity axes may induce a non-uniform convergence 
in series expansions obtained in the original Sinh-Gordon model. Consider for instance 
the total variation of the central charge Ac = c uv — c w going from the short to the large 

2 This is the structure for FF of odd operators. Analogous structure arises for the FF of even operators. 
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distances. For both the SGM and the RTM, Ac = 1. Let us try to express it as a sum-rule 
fulfilled by the two-point function of the trace 0(x) ||, |17| 



Ac = £- J r 2 < e(r)9(0) > d 2 r = £ Ac (2n) , (3.1) 

where Ac( 2n ) is the contribution to the variation of the central charge coming; from the 
2n-intermediate states. In the original SGM with real coupling constant, the convergence 
of the series to the value Ac = 1 is extremely fast and almost saturated by the two- 
particle contribution Ac' 2 -* [11]. This has to be expected, given the massive behaviour of 
the model and the threshold suppression phenomena analyzed in [[H|]. Similar behaviour 
has been also observed in supersymmetric models ||19|| . However, in the RTM the situation 



is drastically different. Consider initially the two-particle contribution to the c-theorem 
sum rule 

The plot of such a quantity (fig. 1) shows that Ac^(/?o) monotonically decreases from the 
value very close to 1 at /3 = (corresponding to the Sinh-Gordon self-dual point) to 1/2 
for (3q — > oo. The asymptotic value 1/2 can be easily obtained analytically by noticing 
that 

E((1,p ) = wah^h(P,Po) , (3-3) 




and therefore for /3 — > oo the integral (|3T2| ) simply reduces to 



, . 3 [°° Ji0 sinh 2 /? 1 , n A . 

Ac^W-oc)^-/ (3.4) 

Concerning the higher particles contributions Ac^ 2n \ all of them vanish in the limit (3q — > 
oo. In fact, the 2n-particle FF entering the formula ( |3.1|) for Ac < - 2n ' is given by eq. (|2.9| ) 
and after the analytic continuation they may be written as 

Ftoi/h, ...,&) = 27rm 2 g 2n ((3 ) Q 2n (l) ]J ^^^'^ , (3.5) 

i<j x * ' ^3 

where g2n(Po) = (4 cosh/? ) ?l ~ 1 A/' 2n ^ n ~ 1 ' ) (/9o)- Analogously to the two-particle case, the 
/^-dependence coming from h(/3ij) is strongly suppressed in the integration over rapidi- 
ties and the asymptotic behaviour in /3q of Ac' 2 "' is only determined by the exponen- 
tial factors contained in g 2n and Q2n(l)- I* 1 the large (3 limit, jV((3 ) ~ exp (— ^-J 
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and then #2n(A)) ~ exp|— (n — l) 2 \(3 \Y On the other hand, for /3 — > oo <22n(l) ~ 
exp {(n — l)(n — 2)|/? |} "P{xj) where V(xi) is a symmetric polynomial. So, for n > 1 
Ac( 2n )(|/3o| — > oo) — > as exp (— (n — l)/5 )- Therefore the result of the series (|3.1| ) is 
Ac =1/2 instead of Ac = 1, i.e. a violation of the c-theorem sum rule. 

Although striking, the non-uniform convergence of the series has a natural interpre- 
tation once the nontrivial interplay between the two scales /3 and (3q of the problem is 
correctly taken into account. In fact, since the n-particle contribution in ( |3.1| ) behaves as 
e -n(Mr) ^ gj ven an y i en gth scale r there is always an integer N r such that the states with a 
number of particles n > N r give a negligible contribution to the series (|3.1|) . This means 

N 

that any partial sum Acn = Ac' 2 " 1 -* only reproduces the variation of the c-function 

m=l 

from the infrared limit r = oo up to a certain scale r^ N '. In usual situations, when c(r) is 
a smooth function in the deep ultraviolet region, the first few A C ( 2n > are sufficient to give 
the correct value of Ac, with high level of precision. But for the RTM this is not the case. 
Consider a scale r\ such that c(ri,/? = 0) > 1/2 (fig. 2). According to the results of the 
TBA analysis, after the first jump from to 1/2, the function c(r, /?o) stays constant at 
1/2 until a value proportional to e - '*'/ 2 is reached and, only at this point the second 
jump takes place. The other jumps occur at r n ~ e -\Po\(n-i)/2 an( j f Qr ^ _^ ^ they 
accumulate to the origin. Truncating the series (|3.1| ) to any N, there is always a value /3q 
such that c(r[ N \ |/3 | > = 1/2; i-e. the point of the first jump is always ahead of the 
corresponding length scale r[ , however small r[ N ^ may be, and therefore 

lim lim Ac N (p ) = \ . (3.6) 



4 Collapse of the Sinh-Gordon Model to the Ising 
Model 

Taking the limit /3 — > oo (keeping (3 fixed), the S- matrix of the SGM goes to S = —1, i.e. 
to the S'-matrix of the thermal perturbed Ising model. Together with ( |3.6|) , these results 
naturally suggest that for (3q — > oo the Hilbert space of the original SGM collapses to 
that of the Ising model, spanned in the local sector only by three independent families of 
fields, those of identity {1}, magnetization {a} and energy {e} operators. It is therefore 
interesting to find the mapping between the operator content of the two models. 

It is easy to see that the elementary field <j)(x) of the SGM is mapped onto the mag- 
netization operator a(x) of the Ising model. In fact, analytically continuing the FF (|2.8| ) 
and taking the limit (3q — > oo, the (3q dependences coming from different terms of the 
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original expression compensate each other and we obtain the following finite result 

2n+l a 

Ft +1 ((3 u . . . } f3 2n+1 ) -> I] tanh^ . (4.1) 

i<j 

These are precisely the FF of the magnetization operator u{x) of the thermal perturbed 
Ising model p0[ ^T]. This field belongs to the interacting sector of the theory and its 



correlation functions satisfy non-trivial differential equations [22, |23[. Notice that in this 
limit the boundary conditions of the field <fi have been modified: in the original SGM its 
FF vanish for large values of $ whereas in the resulting expression ( J4.1| ) they go to a 
constant. 

On the other hand, taking the limit (3q — > oo for the analytic continuation of the FF 
of G Q2.9Q , all of them vanish but F 2 = 2irm 2 sinh/5/2. Hence the operator 0(x) of the 
original SGM is mapped onto the energy operator e(x) of the Ising model. This is a free 
field (a result which is manifest by the absence of higher FF) and its correlators can be 
easily expressed in terms of Bessel functions. Also in this case the boundary condition of 
the field has been changed, since originally F 2 e goes to a constant for large values of /3j 
whereas after taking the limit (3q — ■> oo it diverges at infinity. 

It is also interesting to analyze the behaviour for /3 — > oo of the kernel solutions. In 
this limit the recursive equations (|2.6|) become 

Qn+2(— x, X\, . . . , x n ) = — x n+1 <j n Q n (xi, . . . , x n ) n = odd (4 2) 

Q n +2(-x,x,x 1 , . . . ,x n ) =0 n = even 

The kernel solutions of the Z 2 even operators of the original SGM are therefore mapped 
onto the free sectors of the Ising model, i.e. those given by the identity and energy 
operators. Indeed, their FF are different from zero only at a given level n in the number 
of external particles (where they coincide with S n defined in sec. 2.2.2) and, due to the 
second equation in ( |4.2|) , they decouple from the rest of the recursive chain. Correlators 
of the operators defined by such FF can be also expressed in terms of Bessel functions. 

Such a decoupling in the recursive chain does not occur, on the contrary, for the 
kernel solutions of the odd operators of the original SGM. Their explicit expressions may 
be written as determinants of minors of the matrix E n . In fact, consider the half-infinite 
chain of FF Qi+ 2m ( n °dd and m = 1, 2, . . .) satisfying the first equation in fl4.2p, with 
the initial condition 

QnU = -x n+1 a n Z n . (4.3) 
It is easy to see thatf] Q^+2 = Pn+2] («+!,„._ i) and in general 



(n) 



+2m 



.[£„ 



+ 2mJ( n + 2m-l |n + 2m _ 1 ). 



(2^,n+3). 



(4.4) 



^.n+l) 



We denote by [j4]( a ,&) the determinant of the matrix obtained by A eliminating its a row and b column. 
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Such FF define matrix elements of operators belonging to the magnetization sector. For 
instance defines the FF of the magnetization operator itself whereas those of 
the operator = (o~(x) + 1/M 2 d 2 a(x)) etc. In general such operators have the distin- 
guishing property that their two-point correlation function < O^ n \r)O^ n \0) > decreases 
at infinity as exp[— nMr]. 

5 Conclusions 

The thermal perturbed Ising model is the first model in the staircase series defined by 
the RTM. Using the analytic continuation which links the Sinh-Gordon Model to the 
RTM, we have seen that in the limit f3$ — > oo the elementary field of the SGM becomes 
the magnetization operator a of the Ising model. This is still an interacting field with 
non-trivial form factors. On the other hand, the field and other even operators of the 
SGM are mapped into the free sectors of the Ising model. 

It would be interesting to extend the analysis of this paper to the higher models of the 
RTM and to find the correlation functions of the QFT associated to the corresponding 
massless Renormalization Group flows. 

Acknowledgements. We are grateful to A. Schwimmer for useful discussions. 

References 

[1] Al. B. Zamolodchikov, Nucl. Phys. B342 (1990), 695. 

[2] Al. B. Zamolodchikov, Nucl. Phys. B358 (1990), 497; Nucl. Phys. B358 (1990), 524. 

[3] A.B. Zamolodchikov, Al.B. Zamolodchikov, Ann.Phys. 120 (1979), 253 

[4] A.B. Zamolodchikov, in Advanced Studies in Pure Mathematics 19 (1989), 641; Int. 
J. Mod. Phys. A3 (1988), 743 

[5] G. Mussardo, Phys. Rep. 218 (1992), 215, and references therein. 

[6] Al.B. Zamolodchikov, Resonance Factorized Scattering and Roaming Trajectories, 
ENS-LPS-335, 1991. 

[7] A. A. Belavin, A.M. Polyakov and A.B. Zamolodchikov, Nucl. Phys. B241 (1984), 
333 

[8] D. Friedan, Z. Qiu and S. Shenker, Phys. Rev. Lett. 52 (1984), 1575. 



8 



[9] A.B. Zamolodchikov, JETP Lett. 43 (1986), 730. 
[10] M. Lassig, Nucl. Phys. B380 (1992), 601. 

[11] A.Fring, G.Mussardo and P.Simonetti, Nucl. Phys. B 393, (1993), 413 

[12] A. Koubek and G. Mussardo, On the Operator Content of the Sinh-Gordon Model, 
to appear in Phys. Lett. B. 

[13] A.V. Mikhailov, M.A. Olshanetsky and A.M. Perelomov, Comm. Math. Phys. 79 
(1981), 473 

[14] A.E. Arinshtein, V.A. Fateev and A.B. Zamolodchikov, Phys. Lett. 87B (1979), 389 

[15] F.A. Smirnov, Form Factors in Completely Integrable Models of Quantum Field 
Theory (World Scientific) 1992, and references therein. 

[16] I.G. MacDonald, Symmetric Functions and Hall Polynomials (Clarendon Press, Ox- 
ford, 1979) 

[17] J.L. Cardy, Phys. Rev. Lett. 60 (1988), 2709. 

[18] J.L. Cardy and G. Mussardo, Universal Properties of Self-Avoiding Walks from Two- 
Dimensional Field Theory, UCSBTH-93-12, ISAS-93-75. 

[19] C. Ahn, Thermodynamics and Form Factors for Super symmetric Integrable Field 
Theories, IC/93/144. 

[20] B. Berg, M. Karowski and P. Weisz, Phys. Rev. D19 (1979), 2477. 

[21] J.L. Cardy and G. Mussardo, Nucl. Phys. B340 (1990), 387; V.P Yurov and ALB. 
Zamolodchikov, Int. Jour. Mod. Phys. A6 (1991), 3419. 

[22] B.M. McCoy and T.T. Wu, The two dimensional Ising model, (Harvard Univ. Press, 
1982) and reference therein. 

[23] O. Babelon and D. Bernard, Phys. Lett. B288 (1992), 113. 



9 



